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Abstract 



We study a class of one- matrix models with an action containing nonpolynomial terms. By 
tuning the coupling constants in the action to criticality we obtain that the eigenvalue density 
vanishes as an arbitrary real power at the origin, thus defining a new class of multicritical matrix 
, models. The corresponding microscopic scaling law is given and possible applications to the chiral 

^ • phase transition in QCD are discussed. For generic coupling constants off-criticality we prove 

5^ , that all microscopic correlation functions at the origin of the spectrum remain in the known Bessel 

universality class. An arbitrary number of Dirac mass terms can be included and the corresponding 
massive universality is maintained as well. We also investigate the critical behavior at the edge of 
the spectrum: there, in contrast to the behavior at the origin, we find the same critical exponents 
as derived from matrix models with a polynomial action. 
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1 Introduction 



The study of random matrix models as simplified versions of more complicated interacting Quantum 
field theories has often been very fruitful. For instance, there are zero-dimensional matrix models with 
applications ranging from Quantum Gravity to Quantum Chromodynamics (QCD). For a review on 
these topics we refer to and p|, respectively, as well as to |^] for applications in other fields. In most 
of the cases it has been sufficient to consider the so-called Gaussian matrix model, where the random 
matrices (that describe for example an underlying Hamiltonian) have entries which are Gaussian 
random variables. The reason for such Gaussian matrix models being generic is their property of 
universality. Universality in random matrix models means that when reweighting the matrix elements 
differently within a certain class of perturbations around the Gaussian models the same correlation 
functions of the matrix eigenvalues are retained. Different proofs of universality @, ^ @, |8|, IC , 



11 , using various techniques have been obtained in distinct large- A'^ regimes, where N is the size 
of the matrices. 

In general one has to distinguish between the macroscopic large- limit where the oscillations of 
the correlation functions are smoothed and the microscopic large- limit, where the eigenvalues are 
rescaled such that a particular region of the spectral density is magnified. This region might be the 
origin, the bulk or the endpoint of the spectrum and microscopic universality classes emerge corre- 
spondingly. In fact the universal kernel of orthogonal polynomials that generates all the correlation 
functions is described by Bessel, Sine or Airy functions, respectively. The stability of universality 



under changes of the functional form of the measure has also been tested [l^, 14] and again the mi- 
croscopic and macroscopic regime behave differently. In this paper we will see another extension of 
universality by adding nonpolynomial terms to Gaussian (or more general) weight functions. 

However, this is not the main motivation of the present work. In many applications one has to 
take into account extensions of the Gaussian model and to fine tune in order to reach critical points 
(for a review on phase transitions in matrix models we refer to |l^). There may be several different 
reasons for doing so. For instance, in Quantum Gravity the critical behavior is used to enhance the 
otherwise subleading contributions in the 1/A^^ topological expansion and make them all contribute in 
the double scaling limit (for a review see [|l|). In the study of chiral symmetry breaking in QCD |]l^, ^ 
the spectral density at the origin p(0) is proportional to the order parameter, the chiral condensate. 
Terms of higher order than Gaussian may then be used to achieve a phase transition by requiring 
p(0) =0 1 17]. The corresponding scaling exponents and correlation functions will be those on top of 
the phase transition. Such matrix models have already been investigated in J^] where also universality 
was found within a given multicriticality class. But only very specific transition points could be reached 
in jl7], where p{X = 0) ~ A^™' with m G N. Our aim here is to complete this study and find new 
critical models where the density vanishes with an arbitrary real power. For that purpose we will add 
a nonpolynomial part to a Gaussian potential, or to a more general polynomial potential, and study 
its effect. It has been argued in j|l8| that at the temperature induced chiral phase transition in QCD 
the eigenvalue density precisely behaves as p{X = 0) ~ A^/*^, with mean field exponent 6 = 3. 

Single nonpolynomial potentials have been already studied in J|l^, |2^, |2^. They have been 
shown to fall into the microscopic universality class of the sine kernel in the bulk and at the origin of 
the spectrum. It does therefore not come as a surprise that the extensions of the model we study here 
fall into the corresponding massless and massive Bessel universality classes IC] when Dirac mass 
terms are included. This holds as long as we keep away from criticality or in other words when we 
maintain p(0) 7^ 0. The investigation of new multicritical points is however only possible with several 
interaction terms in the potential, polynomial and nonpolynomial. This is our starting point. 

Nonpolynomial potentials have also been studied in the context of Quantum Gravity 



We will analyze whether the terms we introduce may also change the critical exponents ]P5] in the 



24], 
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appropriate large- iV limit. Finding real critical exponents would correspond to new representations of 
conformal field theories by simple one-matrix models. However, our findings in this respect are negative 
in the sense that the perturbations we consider here will not lead out of the known universality classes. 
Nevertheless it would be very interesting to find realizations of such simple one-matrix models instead 
of more complicated multi-matrix or 0{n) models |26, 27 1. 

The matrix-model approach to the chiral phase transition we described so far is not the only 
possible one. The phase transition can also be studied by adding an external field |18, 28] or by 
allowing more substructure for the matrices such as additional spin or gauge group degrees of freedom 
|29| |. While such extensions are more close to the phenomenology they are certainly more difficult 
to solve exactly without any further assumptions. However, let us emphasize that such a matrix 
model [^] was used to deduce the phase diagram of QCD realistic for two massless quarks in the 
temperature-density plane and that other approaches have confirmed this general picture (see e.g. 
PH). It consists of a second-order line beginning at a critical temperature with zero density that joins 
a first-order line starting from a critical density at zero temperature in a tricritical point. 

In the multicritical matrix models considered in [^] the phase transition is of third order. Therefore 
these models may become interesting precisely at the tricritical point. They may also be relevant for 
the chiral phase transition occurring when the number of flavors Nf is raised above the critical point, 
where the beta-function of QCD changes its sign. The relation of the large-A^y phase transition to 
the one occurring in the Gross- Witten model has been pointed out in |^2|. The extension of to 
nonpolynomial potentials we investigate here may also have a phase transition of different order. The 
fact that the critical exponents of |18| are included in our model suggests a possible second order 
transition. However, the continuity of our results at even integer powers suggests that the transition 
remains of third order. We have not been able to give a definite answer to that question. 

The paper is organized as follows. We will mainly concentrate on Hermitian matrix models also 
called the Unitary Ensemble. This model can be used to study flavor symmetry breaking in three 
dimensional QCD ||3^. In the original approach to chiral symmetry breaking in four-dimensional 
QCD llq] a slightly different matrix model is needed, the chiral Unitary Ensemble. In order to 



translate our results to the chiral ensemble we will simply use the relations derived in |17] and |34] 
between the corresponding orthogonal polynomials and kernels. In Section § we first state our model 
with arbitrary polynomial and nonpolynomial power like terms in the potential. It includes Dirac 
mass terms motivated by the application to QCD. These terms are also studied under the name of 
characteristic polynomials. In the same section we calculate the macroscopic large- spectral density 
away from criticality for potentials containing nonpolynomial parts using different techniques. We also 
prepare the ground for Section^ by recalling the orthogonal polynomial method in the version as being 
reviewed in [35|. In Section^ we present our results about generalized universality. We prove that away 
from criticality the correlation functions in the microscopic scaling limit at the origin fall into the same 



massless and massive universality classes as in ^, 10|, also in the presence of nonpolynomial terms. 
The tool 1 35 1 we use is to derive a universal differential equation for the orthogonal polynomials first 
without massive flavors. Then, in order to include an arbitrary number of Dirac masses we proceed 
along the lines of |12] without any restriction on the number of flavors. Our proof is valid for ensembles 
with chiral and nonchiral unitary invariance and is an alternative to the original proofs ||9|, |l^ for the 
massive ensembles. In Section ^ we come to our main result by tuning the coupling constants of our 
model to criticality. After describing in detail the two simplest examples we introduce a general class 
of new multicritical densities. They behave as p{X = 0) ~ A*^"^ with real noninteger k > 1. For 
2m — 1 < K < 2m + 1, m G N, we also give the corresponding minimal critical potentials containing m 
polynomial terms and one nonpolynomial term. In Section ^ we investigate possible critical points in 
the leading order of the free energy at large- for the simplest of such multicritical models. We find 
a nonanalytic behavior at large- when the spectral density is tuned to develop additional zeros at 
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the edge of the spectrum. The critical exponent '^str = — ^ we find for the free energy is the same as 
known from matrix models with polynomial potentials |25]. Section ^ contains our conclusions and in 
the Appendix we collect some results on the so-called two-cut spectral density, where the eigenvalues 
lie on two disjoint intervals that merge at the transition. 



2 The model and its large- solution 

We start by describing the most general matrix model we wish to work with. In order to make sense 
of nonpolynomial terms we give the model right away in terms of the eigenvalues of a Hermitian 
matrix: 



.00 N ( Nj \ 

/ widXr n (^^' + -/) 1^(^)1' ' (2-1) 

■^-~i=l\ /=1 / 



V{\) ^ Vpoi{X) + Vnonpoi{\) ^Y.Y'^^' + Y.H^\''' ' (2.2) 

j=l ^ i=l 

where 1 < ki < . . . < k^^^ € M\N. Here we have restricted ourselves to 1 < Kj as we do not only need 
a normalizable spectral density p(A) enforcing < Kj, but also p(0) < cxd which requires 1 < Kj, as we 
will see. The Vandermonde determinant A (A) = nfc>/('^fc — A/) originates from the diagonalization 
of the Hermitian matrix. The potential V{\) is chosen to be symmetric for simplicity and in view 
of applications to QCD (as we will show, odd powers are not necessary for obtaining multicritical 
potentials). It has been split into a polynomial part, Vpoii^)-, and a nonpolynomial part, Vnonpoli^)- In 
the latter we have allowed for an arbitrary number dnp of real positive powers Kj. In order to achieve 
a multicritical behavior at the origin we will actually need only one of such terms. However, for the 
large- solution presented in this section as well as for the issue of universality in the next section 
any number of such terms may be present. 

The Dirac mass terms in front of the exponential are added in order to study flavor symmetry 
breaking in three-dimensional QCD (see ||3^ for details). The 2Nf Dirac masses occur in A^^ pairs 
±imf. Furthermore, we have included a massless Dirac determinant raised to some arbitrary real 
power 2a > —1. We will need this term later in Section ^ when we switch from the Unitary Ensemble 



of Hermitian matrices in eq. (2.1) to the chiral Unitary Ensemble which is needed in applications to 
four-dimensional QCD. 

In the following we will determine the spectral density p{X) in the macroscopic large- A^ limit. A 
first argument how to calculate it in the presence of nonpolynomial terms comes from a saddle-point 
analysis. To begin with we observe that the preexponential factors containing the mass terms in eq. 
(|]|) can be formally included into the potential by shifting V{\) V{X) - l/A^ln[|Ap° H^iil^^ + 
m^)]. Hence they are subleading in the large- A^ limit and we can set a = Nj = in the rest of this 
section. Only in the microscopic limit they will become important while the macroscopic spectral 
density does not depend on a,Nf ^ 0. 

The spectral density for a single nonpolynomial potential V^(A) = jAj*^ has been calculated already 



in ||19|, mil HI- In order to determine the full density p(A) for the potential eq. ( P^ ) we make use 
of the fact that the saddle-point equation for the spectral density is linear in the potential. We can 
therefore obtain solutions of the general problem by simply adding up solutions for the pure polynomial 
case with solutions for the nonpolynomial case. In order to see that in more details we repeat some 
well known facts about the saddle-point solution of matrix models (for a review see e.g. [Q]). Defining 
the free energy as J- = —ln[Z]/N'^ and switching from discrete eigenvalues to a continuous density 
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one has 

T = I d\p{X)V{X)- I dXdpp{X)p{iJL)\n\X- + K [ / dXp{X)-a] . (2.3) 



We have introduced a Lagrange multiplier K to enforce the normalization of the spectral density p{X) 
to a, which may differ from unity. Furthermore, we consider here the case where the support a of 
p{X) is a compact interval, the so-called one-cut solution. For the two-cut solution we refer to the 
Appendix ^ The saddle-point equation then reads 

= V{X) -2 I dp.p{p)\n\X- p\ + K , AGo- = [-c,c], (2.4) 

J a 

which is often written also after taking the derivative w.r.t. A as 

.„i« = m, .e.. (2.5) 

A — fl Z 

and ^is the Principal Value of the integral. The Lagrange multiplier K can be determined by choosing 



a particular value of A in eq. (2^), e.g. A = S o". Because of the linearity of eq. ( p^ ) two positive 
solutions pi and p2 for two given potentials Vi and V2, respectively, can be simply added under the 
assumption that they have the same support a. This gives a new solution p = pi+ P2 for the potential 
V = V\ + V2 which is nonnegative on the same support a. Thus we can take known solutions of the 
saddle-point equation for polynomial and nonpolynomial potentials, fix their support to be equal and 
add them. The normalization to unity of the new spectral density, p{X) = PpoiiX) + PnonpoiWi has 
to be imposed eventually. This condition uniquely determines the dependence of the end-point c on 
the coupling constants in the potential or, which is the same, it gives a relation among the coupling 
constants for any given fixed positive number c. 

Let us give an example. Choosing the simplest polynomial potential, i.e. a Gaussian potential 
VpoiW = ^52 A^, the corresponding spectral density is the "Wigner semi-circle" 

PGaussW = |^\/c2-A2 , A G [-C, c] . (2.6) 

If one fixes the normalization condition then (72 and c must be related to each other. In particular, the 
normalization condition of the density to unity reads \g2(? = 1- Here we consider c as an arbitrary 
given positive parameter. For the nonpolynomial part we choose VnonpoiW = 9|A|'*, which is called 
Freud weight in the mathematical literature. The corresponding spectral density normalized to unity 
has been calculated in [^, ^] (see also references therein). The un- normalized density sharing the 
same one-cut support [— c, c] will be re-obtained in this section and can be written as (see also [21]): 

PFreudW = ^'^^'^ / , ds ^ (2.7) 



TT V 2' ' 2' C2 / V C2 



gKHc'^-^^f^ K 3 ^ A2\ / A2 ^ ^ 

-F — U/l-— , Ag[-c,c], (2.^ 



where 



r(f 



1- 



B is the Beta function and F is the Hyper geometric function. It can be easily seen from the integral 
that for < K the density is normalizable and that for 1 < k it holds p(0) < oo. For a discussion 
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of < K < 1 we refer to 21 1. Using some identities among Hypergeometric functions for k ^ 
1, 3, 5, . . ., eq. (|2.8|) can also be written in the useful form: 



PFreudW 



gnHd 



K-l 



/ 1 


' KTT' 




tan 




(^) 


I 2H 


I 2 \ 






(2.10) 



Imposing the condition 1 = Hgc'^ we would obtain a normalized density, reproducing |19, 22 1. By 
adding the two solutions eqs. ( |2.6D and (^]^) which now have the same support, we obtain for the full 
spectral density 



P(A) 



PGauss (A) + PFreud (A) 



2tt 



IT 



ds 



A2 



A e hc,c] 



(2.11) 



By imposing dXp{X) = 1 we obtain the following relation between the coupling constants and the 
endpoint c of the support: 



1 



-020^ + Hgc" 



(2.12) 



This relation can be easily checked by comparing to the known polynomial potential with quadratic 
and quartic coupling [36| by setting k = 4. We can immediately draw two consequences. First, both 
parts of the density, pcaussW and PFreudW have an expansion in powers of around the origin 
with the latter having an extra term ~ A''^^ as one can see from eq. ( |2.1C ). Consequently a critical 



behavior at the origin can be achieved by tuning the coupling constants §2 and g such that the leading 
constant term, /o(0), vanishes. This will be analyzed in great detail in Section ^. Second, it is easy 
to see how to generalize the above adding of spectral densities. Since the density for a polynomial 
potential is already know to most generality [Q] we can use the above procedure to add successively 
single nonpolynomial potentials of higher and higher real power to obtain the density for the potential 
eq. (^. 

A standard way to obtain eqs. (|2.7D and (|2.8|) , is by introducing the resolvent G{z): 



G{z) 



z — n 



z e C\a 



(2.13) 



in terms of which the saddle point equation ( ^ ) reads G(A + iO) + G(A — iO) = ^'(A), X G a. Clearly, 
G{z) has a cut at z £ a. The (unique) solution which is bounded at the end points and which behaves 
like G{z) ~ 1/z at large \z\ is (see e.g. 



G{z) 



dt 



V'it) 



27r {z-t)V^ 



1 + K 



2 ■ z2 y V C2 



(2.14) 



in the one-cut case, for V{X) = ^lAj". For any real positive k, the spectral density eq. ( |2.8| ) is then 
recovered by means of 2-Kip{X) = G{X — iO) — G{X + iO), X G [— c, c] (after applying an identity for 

2 2 

the Hypergeometric function corresponding to the transformation % — > 1 — The normalization 



condition is readily obtained from the large-2: asymptotic behavior of eq. (2.14). In fact by using that 
-F(a,/3,7,0) = 1 one has G{z) ~ gHc'^/z = 1/z which is the previously stated formula gHc'^ = 1. 
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In the remaining part of this section we obtain the spectral density by using another technique, the 
method of orthogonal polynomials. It will turn out to be very useful also in the next section where we 
prove the universality of all microscopic correlation functions away from criticality. Let us begin by 
repeating a few well known formulas [^] in order to fix our notation. We choose a set of orthonormal 
polynomials, 

r dAz/;(A)Pi'^^ '")(A) Pi''^^ '")(A) = , (2.15) 

J — oo 

where the weight function is taken from eq. (2.1): 



Nf 

w{X) ^ |A|2" J] (a2 + mj) e-^^W . (2.16) 
/=i 



The fe-point correlation function of the eigenvalues, 

N N 



i?(J^^'")(Ai,...,Afc) = ^^^^ 1^ n dX, l[w{X,) |A(A)|2 , (2.17) 



i=k+i j=i 



defined as the {N — A:)-fold integral over the integrand of eq. ( |2.lD , can be expressed through the 
kernel of the orthogonal polynomials p8[ 



i?r^--)(A„...,A.) = (^,Jet^J^r^"^(A.A,)] . (2.18) 



The kernel is given by 

N-l 



K^'-hx,,) = ^.(A)i.(,)i ^pf^^-)(A)/^(^^^-)(,) 



1=0 

{2Nf . {2Nf ,a) , , (2Nf ,a) , , {2Nf ,a) , . 

= w{X)2w[ri)'2 Tat , (2.19) 

A — ?7 

where in the last step we have used the Christoffel-Darboux identity. It follows from a three step 
recursion relation, which any set of orthogonal polynomials obeys. For any symmetric potential V{X), 
polynomial or not, it reads 

APr^-)(A) = r^^Afr'W + ry:T\x) , (2.20) 

under the assumption that the integrals over the P^^^'°'\X) exists. The recursion coefficients are 
determined from the so-called string equation 

n = -r„ r dAt.'(A)Pi'^^'")(A)pf_^^'")(A) . (2.21) 



In the case of a = Nf = 0, which is all we need in order to determine the macroscopic spectral density, 
one can show that by setting = 2c eq. ( |2.21 ) is equivalent to the normalization condition of the 
spectral density arising from the saddle-point analysis (e.g. see eq. ( 2. 12] )). 

In order to determine the spectral density one has to consider also differentiation of orthogonal 
polynomials. Here we will use the Shohat method, which is reviewed in great detail in |^^. It 
consists in a particularly useful rewriting of P^{X), which is in general a linear combination of all 
lower polynomials, in terms of two auxiliary functions An{X) and Bn{X). 
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For the rest of this section we will restrict ourselves to the massless case, Nf = 0, (and suppress 
the index) . The massive case with A'^^ 7^ will be treated later in Section |3|. Furthermore we absorb 
the massless prefactor |Ap° from the weight function into the potential 



K(A) = F(A)-^ln|A|. 



Defining the auxiliary functions[| 

An{X) = Nr, 



dt w{t) 



t2- A2 



(2.22) 

(2.23) 
(2.24) 



it is easy to show using the Christoffel-Darboux identity that the following relation for differentiation 
holds for any finite n and for any potential Ki(A) 



'^-^^ = ^n(A)PS(A) - i?„(A)p(-)(A). 



Furthermore these functions obey the following identity at any finite n [35|: 

A 



i?„(A)+i?„_i(A) + iVy^(A) 



rn-l 



-An-l{X) . 



(2.25) 



(2.26) 



The relation ( 2.25| ) will be the key to derive an exact differential equation for the wave functions 

4°)(A) ^ «;(A)lpW(A) (2.27) 

in the next section and to prove universality. 

Now we are using eq. (|2.25| ) to derive the spectral density for the general potential defined in eq. 
(p.2|). From eqs. ( 2.18| ) and ( p. 19 ) the one-point function or spectral density is simply given by a 
single kernel at equal arguments: 

(2.28) 



= w{\)r^{pP\x)PP_,{\)-PP\\)PP^,{\) 
Using the recursion relation ( |2.20|) and the identities ( 2.25 ) and (|2.26| ) we arrive at 

^7v(A)^{,li(A)2 + ^An-i{X)^^S\^? 
rN-i 

(^-±-An-,{X)+Bn{X) - S;v-i(A)) V'S'^(A)vSv"li(A) 



<(A) 



(2.29) 



which is exact for any finite A^. We will now take the macroscopic large- limit where we smooth the 
oscillations of the correlators. By considering smoothed moments of the orthogonal polynomials at 
large- the following result has been shown in [^]. For any set of orthogonal polynomials of a given 
measure w(A) which first leads to a large- A^ spectral density with single interval support a and second 
has recursion coefficients obeying the condition tat-i-i -1-2.... — > r^r = 2c when A^ — > 00, it holds 

#(A)2 = l^^9{c'-X'), 



vr Vc2 - A2 



C^(A)^]v!i(A) 



1 



A 



cvr \/c2 - A2 



A^). 



(2.30) 



^Compared to |33] we have already used the symmetry of the potential (see also appendix A there). 
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Here the bar denotes the smoothed large- A'" hmit. Let us stress that the derivation of eq. ( 2.30I ) in the 
Appendix A of is purely algebraic. No assumption about the analytic structure of the spectral 
density in the complex plane has been made, in particular no square-root behavior has been assumed. 
Using these results we can take the same smoothed large- limit for the functions A]\;(X) and B]\i{X) 
defined in eqs. ( p. 23 ) and ( 2.24 ). Smoothing inside the integrand we obtain 



^(A) ^ -An{X) 



B{X) ^ -Bn{X) 



c r tv'jt) - xv'jx) 

Jo - A2 

1 XV'{t)-tV'{X) 

IT 



1 



2_-f-2 ' 

t 



t2 _ A2 V^^3t2 



(2.31) 
(2.32) 



Note that compared to eqs. ( |2.23| ) and ( |2.24| ) we have dropped the index a in the potential V{X) or 
in other words the term 2a In |A|/A^ from eq. ( |2.22 ). This can be seen as follows. Because of parity it 
drops out in ^Ar(A) already at finite N. In one can split Bn{X) = B^^ reg(A)-|-(l — (— l)^)a/A 
where the first term only contains V instead of Va- The latter term obviously has no smooth limit 



but it is suppressed by in eq. (|2.32[) . 

We can now apply the smooth limit to the density eq. ( p. 29 ) itself and we obtain 



(A) 



— A(A)v'c2- A2 , A G 



-c, c 



cvr 



(2.33) 



In particular it holds that 



1 



p(0) = - A(0) 



vr 



(2.34) 



Eq. (2.33) together with eq. (2.31) constitute our main result of this section. Let us stress that it 
does not necessarily imply that p{X) has a simple square-root cut as a function in the complex plane. 
The function ^(A) having nonpolynomial parts inside the potential may precisely cancel the factor 
Vc^ — A^, as we shall see in a simple example now. Let us choose again the nonpolynomial potential 



Freud 



(A)=5|Ar 



1 < K G 



(2.35) 



From eq. ( 2.31 ) we obtain for the auxiliary function 

t'" - A" 1 



A{X) 



c 
vr 



dt 







XWc 



2\/? 



A2 



.1 K K, C 

kA2 I ' 2 + 2' ^ 2' A2 



(2.36) 



In order to be able to expand the density around the origin and to see the square root cancelling in eq. 
( p. 33 ) we still have to use some identities for Hypergeometric function (and to analytically continue 

to A2 < c2) 



1/1 K K c2 \ 

A2 r' 2 + 2' ^+ 2' A2 j 



1 1 1 K 

TW^X V2'2'^+2'A2j 

1 /^r (i + f)r(i-f) ^,.,A 



K-l 



VA2 - c2 



/vr 



r(i + f)r(f -1 



r(f)r(f + i) 



c 12' 2 



K 3 

2' 2 



K X' 



(2.37) 
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Inserting this back into eq. ( 2.36| ) we obtain with eq. ( 2.33 ) the final expression for our density as 
given in eq. ( 2.10| ): 



PFreud{X) 



gtiHc 



vr 



2H 



tan 



K-1 



+ 



-F 



1 1 

2' 2 



K 

2' 



(2.38) 



This solution for the density we have derived can be seen to be equivalent to the first representation 
given in eq. (^^), by using the identity relating Hypergeometric functions with argument z and 1 — z 
respectively^. The advantage of eq. ( 2.1C| ) is that we can immediately read off the expansion of the 
density at the origin. For noninteger k > 1 it consists of a single term ~ A'^"^ plus a power series 
in A^. It is the first term which will be responsible for the new critical behavior to be discussed in 
Section |^. 



3 Microscopic universality 



In this section we prove that in the microscopic large-A^ scaling limit at the origin all correlation 
functions of our model ( |2.lD are universal, as long as the couplings in the potential V{\) eq. (2.2) are 
generic (=noncritical) and we stay in the phase of a single interval support. They belong to different 
universality classes depending on the number of massless and massive flavors, a and Nf respectively, as 
well as on the rescaled masses ruf. The universal parameter will be /9(0) which encodes the influence of 



all the coupling constants in the potential via eq. (2.34). Here we will make use of the powerful results 
on universality for polynomial potentials which have been obtained for massless ^, ^, ^] and massive 
|9|, |lO|, ^ unitary and chiral unitary matrix models. Our extension to nonpolynomial potentials works 
as follows. We start with the unitary model without massive flavors, Nf = 0, resuming the method 
of orthogonal polynomials used in the last section. A differential equation for the asymptotic of the 
polynomials will be derived for arbitrary a > — 1 following B. Once we have established massless 



universality we can use the method of |12] to extend it to an arbitrary number of massive flavors A'^^. 
In contrast to the chiral case no restriction on Nf to be even is needed here, as has been noted already 
in the extension to the non Hermitian case [40|. In a final step we use a direct relation between the 
polynomials of the unitary and the chiral unitary ensemble to prove universality also in the chiral 
case, for arbitrary a and Nf. 

Let us start with the unitary case and Nf = as in the end of the last section. Following |35|, eq. 
( p. 25 ) for differentiation on the polynomials can be used together with the recursion relation ( |23oD 
and the identity ( 2.26 ) to derive an exact second order differential equation for the wave functions 



i^ri^iX) eq. ( p. 27 ) for any finite n (for the present form see ||T 



V^i-)(A)" - F„(A)V^i")(A)' + Gn{\)M^) 



(3.1) 



where 



Fn{\) 
Gn{\) 



<(A) 



An{\) ' 



rn-i 
+ 



An{\)An-l{\) 



N , 



Bn{\) + yF^(A) 



Bn{\) + yKi(A) 



A„(A) 



N ^ 



Bn{X) + -V:,{X) 



(3.2) 



When K is not an odd integer number. For k an odd integer number the formula is shghtly more complicated |g9|. 
However, formula (2.1C) is completely fine here because as we will see in Section^, odd integers k do not affect the issue 
of multicriticality we are considering. 
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Imposing a smooth limit for the recursion coefficients tat as well as for A]\f(\) and the regular part of 
Biy{X), the function can be entirely expressed in terms of the limiting function ^(A), due to 

the identity ( ^2^ . We obtain 



(3.3) 

The microscopic large- limit at the origin is defined^ as 

( = NX, (3.4) 

where ^ is kept fixed at large- A^. This rescaling is the only consistent way to obtain a finite large- 
differential equation from eq. (|3.3| ) and it agrees with the rescaling in the pure Preud case obtained 
in |19, 2C, 21, 22| when translating to their conventions. Eq. ( |3.4D holds despite the fact that ^'(A) 



diverges ~ A" ^ for 1 < k < 2. Even in that case the logarithmic derivative of A{X) (as well as 
A^yl(A)^) is suppressed and we arrive at 



i^NiO" + (^(0)' + ^^^2 J ^NiO = . (3.5) 

This equation of Bessel type entirely determines the asymptotic large- behavior of the wave function, 
distinguishing between even and odd polynomials. Moreover it is universal since it only depends on 
the potential V{X) eq. (2^) through the universal parameter ^(0) = TTp{0)- We refer to ||6|, |8| for 
the explicit solution of the differential equation as we are only interested in proving an extension of 
universality to nonpolynomial potentials here. Similarly we refer to Q for the correlation functions, 
where detailed expressions can be found. We only wish to stress here that from eq. (|3.5| ) the rescaled 
microscopic kernel 

= hm -K^'^U^ ^ 



/^r(e,C) - lim ^i^ri^,^) (3.6) 



follows as it can be entirely expressed in terms of the asymptotic wave functions (see eq. ( p.l9D ). The 
microscopic correlation functions defined as 

can then be obtained from the rescaled version of eq. ( |2.18| ). 

Up to now we have determined all microscopic correlation functions from eq. (|3.5| ) for an arbitrary 



a > —1 and proved their universality. In the next step we reintroduce the mass terms from eq. (2.1) 
and determine their universal correlations. Following [^] these masses can be put into the Vander- 
monde determinant leading to a relation between the massive A;-point correlator and the massless 
(k + A'^^)-point correlator analytically continuated in some of the arguments. 

To be more precise we can rewrite the mass terms times the Vandermonde in the integrand of eq. 
(0) as 

nn(A,-.»,)(A, + .»,)|A„(A„...,A„)P = |A«.MAl.....A..imi... 

IAtv (zmi,...,«miv )|2 



^For 1 < K < 2 the slope of p(A) is infinite at tlie origin (see Figure ^) and thus the eigenvalues have to be unfolded 
before rescaling. We thank K. Splittorff for raising this point. 
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Here we have explicitly given all the arguments of the Vandermonde determinant. Because of the 
pairing of the 2Nf masses in pairs ±imf coming from QCD3 we do not need to impose an extra 
degeneracy to the power of the Dyson index /3 = 2 as it was done in the chiral case [||]. This 
observation valid for the unitary case was made in [40|. Inserting eq. ( |3.8D into the definition (2.17) 
we see that the right hand side is proportional to a massless k + A'^y-point function, with no mass 
terms. Inserting all normalization factors we arrive at |40|^ 

(Ai,...,Afc) = (o,„) ,. : ^ • (3-9) 

Since for the right hand side we have already shown its universality in the microscopic limit the same 
statement holds for the left hand side. Here the masses have to be rescaled the same way as the 
eigenvalues, keeping /x/ = Nnif fixed. We have thus calculated the massive correlation functions and 
at the same time given a proof for their universality. Explicit expressions in different equivalent forms 
can be found in |1C, 41 1. We have thus not only extended the universality proof of |l^ to potentials 
of the form in eq. (2^) but also considerably shortcut proof and calculations. 

After having completed the universality proof for the full model eq. ( |2.1| ) it remains to repeat the 
analysis for the corresponding chiral model. The partition function can be written in two ways: 

({-/}) = Tn (d^^ (^^ + -/) e-^^— |A(A)P 

+ mj) e-^^^'^-'(^?)j |A(A2)|2 , (3.10) 

VcMral{\^) ^ ^^A2^- + ^2/i,|Ar^ , (3.11) 

i=i j=i 

where in the second line we have simply changed variables A ^ A^ . In order to calculate the correlation 
functions we introduce again orthogonal polynomials in analogy to eq. (2.15): 

/"CXD J 




By comparing to the unitary orthogonal polynomials eq. ( p. 15 ) we can immediately read off that the 
two sets are related. This relation has been already found in [|l^] and trivially extends to the massive 



case 



3|]: 

^2N - ^N,chiraiy^ I ' {6.16) 

where we also have to identify the corresponding potentials 

2V{\) = VcMral{\^) ■ (3.14) 



*To be precise the relations holds for finite- TV only in the normalization of At large- A'^ this becomes irrelevant. 
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In other words the even subset of the polynomials of the Unitary Ensemble at shifted a + ^ is sufficient 
to construct the full set polynomials of the chiral ensemble. The same is true for the wave functions. 
Since we have started with a real a this shift is possible. We can thus borrow the full machinery for 
the asymptotic of the polynomials in the unitary case, including universality. We have shown already 



that through eq. (3^) the asymptotic of the wave functions in the massless unitary case are the same 
for polynomial and nonpolynomial potentials and that they are universal. Through eq. ( |3.13| ) the 
same statements immediately translate to the wave functions of the massless chiral ensemble and thus 
through eqs. ( 2.19| ) and (2.18) to all massless correlation functions. For the massive correlations we 
have to invoke another relation since for the Unitary Ensemble we did not determine the massive 
orthogonal polynomials. We directly determined the massive correlation functions through eq. (|3.9D 
and thus we cannot use eq. (3.13) to read off the massive chiral orthogonal polynomials. Fortunately 
the corresponding massive kernels of the two ensembles, chiral and non chiral, can be directly related 
as well by using eq. ( 3.13| ) p^]: 



K 



{Nf,a) 



- K. 



{2Nf,a+l/2) 



2N 



(A, 7?) + K, 



(2Afj,Q+l/2) 



2N 



-A,r?) 



(3.15) 



The massive kernel of the Unitary Ensemble on the right hand side follows from the massive two-point 
function that we have already determined in eq. (|3.9| ). Inspecting eq. ( 2.18 ) we see that it follows 
from 

1 



K 



{2Nf,a) 



N 



-R 



{2Nf,a) 



'N,con' 



where the connected part of the correlation function is defined as 



{2Nf,a) 
N,conn 



(A,r?) 



R 



{2Nf,c 



N 



'(A,r?)-i?: 



{2Nf,a) 



N 



(A)i?r^-^ 



(3.16) 



(3.17) 



Using again eq. ( 2.1g| ) we have finally also determined all massive correlation functions of the chiral 
ensemble. Our consideration establishes an alternative proof to the original one in |p, in a considerably 
simpler way. This is the gain we made by first solving the nonchiral ensemble. Let us stress that in 
contrast to [|l2| we have no restriction on the flavor number A'^^. 

Finally we mention that the identification of the potentials, eq. ( 0.14 ) is particularly useful when 
studying the multicritical point as it will be done in the next section. We can therefore restrict 
ourselves again to the unitary case. 



4 New multicritical behavior at the origin 

Multicritical points occur when the spectral density p{\) develops additional zeros inside (or at the 
edge of) the support of the eigenvalues. Since in this section we are interested in applications to the 
chiral phase transition in QCD, with the chiral condensate being proportional to /o(0), we restrict 
ourselves to additional zeros at the origin of the spectrum. The multicritical points can also be 
thought of the support merging together from two (or more) segments. In fact for the question of the 
order of the phase transition the free energy and its derivatives have to be calculated from both sides 
and compared. In general these transitions are of third order which also holds for additional zeros 
developing elsewhere on the support. 

Before defining the new class of critical potentials let us briefly review what is known about 
multicritical models with polynomial potentials. Since we only study the behavior of the macroscopic 
density for the degree of criticality we can set again Nj = a = since it does not depend on these 



parameters. From eq. (|2.3l| ) we deduce that for polynomial potentials the density (|2.33[) consists 



of an even polynomial times a square root ^/c^ — Expanding p{X) at the origin and tuning the 
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couplings for p{0) to vanish we can only achieve that p(A = 0) ~ A^™" for m € N. For the mth 
multicriticality we have to tune precisely m coupling constants. In a set of minimal potentials 
was given. The corresponding mth multicritical correlation functions were determined in and 
shown to be universal. We will find a similar situation here (although we will not touch the issue 
of universality), namely for a spectral density vanishing like a real power p{X = 0) ~ A'"-\ with 
2m — 2 < K — 1 < 2m we have to tune m coupling constants. 

To see this let us start with the simplest example, a Gaussian plus a single real positive power: 



V{\) = + g\Xr . 



(4.1) 



Instead of the solution eq. ( 2.11| ) we use the equivalent representation for the integral, eq. ( 2.1C| ): 



P(A) = 1^ tan 



1) 



2' 2 2' 2 



^:-,^) + ^92^0^, (4.2) 



2tt^ 



where we now have to restrict to 1 < k 7^ 3, 5, 7, ... . The normalization condition eq. ( |2.12| ) which 
we repeat here for completeness fixes one of the coupling constants 



1 



1 



g2C^ + Hgc'' . 



(4.3) 



When expanding the density at the origin the Hypergeometric function as well as the square root both 
have an expansion in A^ and we obtain 



P(A) 



gK 
27r 



tan 



KIT 



A'^-i + 



1 



2k 



gHc''-^ + 52C ) + 0(A 



27r V K - 1 



By tuning the remaining coupling constant we can achieve p{0) = through 



2k 

K-1 



However, only for /t < 3 the first term in eq. (E^) with a real power is dominant and we have 



gK 

0) ~ — tan 
' 2lT 



KTT 



L 2 



1 K — 1 



+ 0{X^) for 1< K < 3 



(4.4) 



(4.5) 



(4.6) 



Consequently the potential eq. ([4. ID together with the condition eq. ( [4.5[ ) and the normalization eq. 
( [4.3| ) defines our first new class (|4.6| ) of multicritical models. A few comments are in order. First of all 
it is easy to see that equations ( [4.5[ ) and ( |4.3| ) have a unique solution with 52 > 0, < for 1 < k < 2 
and with 52 < 0, 5 > for 2 < k < 3. In both cases the slope of p{X) in eq. (4.f:) is positive (infinite 
or finite), as it should be. An example for each of such cases is given in Figure^. 

For real k > 3 the quadratic term in (4.4) is leading and we are back in the class p{X = 0) ~ A^™ 
for m = 1 [|l^]. If we add a higher order term to the polynomial part of the potential, like a quartic 
power for example, we can still achieve that p(0) = 0. There is one free parameter left and we remain 
in the multicritical class of [0] unless we also cancel the C(A^) term. In the latter case, which we will 
give as a second example below we may arrive at higher criticality, p{X = 0) ~ A""^, with 3 < k < 5. 
Before doing that let us remark that instead of eq. (4J) we could have started with a potential 
V{X) = hi\X\'^'^ + /i2|A|''^ instead of perturbing around the Gaussian. Going through the derivation 
again, where we now add two densities of the type eq. ( |2.10| ), we could again achieve that p(0) = 0. 
The leading order term of p{X) would then be to the power of min{Ki — 1, k;2 — 1) 2}. Hence in case of 
1 < Kj < 3 for at least one of the Ki we would be again in class ([4.6D, otherwise in class m = 1 of IT! 



Adding higher order terms, polynomial or not, would then allow to obtain higher criticality as well. 
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-1 -0.5 0.5 l_x -1 -0.5 0.5 1_\ 

Figure 1: The multicritical density for 1 < k = | < 2 (left) and for 2 < k = \/5 < 3 (right). 



In order to keep the algebra simple we make the minimal choice of one real power plus a polynomial 
part for the minimal critical potentials. As an example for the second multicritical class we thus add 
a quartic term to eq. ( [4.1[) and start with 



V{X) = \g2\^ + \gA' + g\\r 



(4.7) 



The density can be obtained using eqs. ( |2.33|) a nd ( 2.31| ) or by simply adding the well known density 
of the purely quartic potential [^ ] to eq. (|2.li| ) and normalizing properly. We obtain 



pW = 1^ tan 



I K— 1 



11 K 3 



gnHc ^ . 
■k{k-1) V2' 2 2' 2 



A2 



2' 



with the normalization condition 



1 = 752c2 ^ 3 ^^^4 ^ ^^^^ ^ 
4 lb 

Expanding at A = and setting the terms of 0{1) and 0{}?) to zero we obtain 



L 2 



A 



with conditions 








2k 

At- 1 

K 



+ ©(A^) for 3 < K < 5 



1 



(4.8) 



(4.9) 



(4.10) 



gHc'' + g2c' + -54C* , 



1 3 
-(7i?c'' - -g2C^ + -gic^ 



(4.11) 



After giving two explicit examples we now wish to come to the general mth multicritical density 
and potential. In order to have the first term in eq. ( |4.8D dominating for 2m — 2 < k < 2m we have 
to add a polynomial potential with at least m coupling constants in order to cancel the first m terms 
in the expansion of the Hyper geometric function there. Generalizing the cases eq. ( [4.1D and eq. ( ^.7| ), 
hereafter we consider the potential 



v{X) =yp„;(A)+5iAr = x^f^A^^- + 5|Ar 



(4.12) 
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The general one-cut solution for the polynomial part is know to be 



PpolW 



7=1 k=0 ^ ^ 



which we rewrite as 



with 



m—l 



PpolW = ^aj-^Jl 



j=0 



A2 



-1 



m— 1 

k=j 



2{k-j) 
k-j 



(- 



2k 



J = 0, . . . , m - 1 



(4.13) 



(4.14) 



(4.15) 



The full density is obtained by adding the Freud-type density eq. (|2.10| ) to eq. (|1|), and the 
normalization condition reads 



m—l 



j=0 



(2i - 1)!! 
2J+i(j + l)! 



(4.16) 



with (—1)!! = 1. When expanding the Hyper geometric function at A = we further obtain m equations 
by setting the terms 0(1), . . . , ©(A^™"^) to zero. We find| 



-1 3 



n 



{21 - k) 



7r(l - k) (2/ + 1 - k) 
The mth multicritical density is thus reading 



j = 0, . . . , m - 1 . 



(4.17) 



P(A) 



gnHc 



7r(«;- 1) 

gn 
+ — tan 
27r 



11 K 3 K A" 

2' 2 ~ 2' 2 ~ 2' 



'^^^J2^-^ A^ 
^^\{21 + 1-k) 



KTT 



\ K— 1 



j=0 I 

2m - 1 < K < 2m + I , 



A2 

^2 



(4.18) 



and it behaves like p{X = 0) ~ A""^ + ©(A^™). The critical values of the m + 1 coupling constants 
g2j and g^ in the potential eq. ( [4. 121) can be obtained by iteratively solving the linear set of m + 1 
equations given by eq. (4.15) and eq. ([4.16 ), which is already of triangular form. Actually there exists 
a much more compact form for eq. ( [4.181) . It can be obtained by iteratively using the identity 



(^-2) 2 p(l 3_K _ 5 
(k-3) ^ V2' 2 2' 2 



2' 



F -, ; ; 

' 2' 2 2' 2 2' 



m 



(4.19) 



1, in 



which can be derived using standard formulas [39| . Going back to the first criticality 
eq. (4.2) the square root and the Hypergeometric function can be put together using the identity 
( [4.19 ) and eq. ( [4.5[ ) (see also eq. ( A.1C| )). At m = 2 for eq. ( [4.8D we can do the same. The resulting 
Hypergeometric function can be simplified again together with the term A^-y/c^ — A^, using the identity 
( [4.19 ) at K ^ K — 2. Proceeding by induction with the general form eq. ( [4.18[ ) we arrive at our final 
result: 



p(A) = — tan 



KIT 



gKHc'^-^ A^"^ 



K 3 

2' 2 



1 ^ 
2' c- 



{21 - k) 



1=1 



{21 + 1- k) 



^The meaningless product HiLi which appears when j = 0, is set equal to 1. 
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Here the vanishing at A = is manifest. Notice that in our parameterization, the end-cut point c is a 
given positive real number: that means that in general one has at most a critical line in the space of 
the parameters of the model. This finishes our complete classification of new multicritical models. 



We turn now to the microscopic scaling limit of the multicritical models. While the differential 
equation ( 0.1| ) together with (|3.2| ) for the wave functions ip^^X) still holds exactly for finite the 
large-A^ limit will change. We first determine the appropriate microscopic scaling limit that corre- 
sponds to eq. (|3.4D at criticality. For that purpose we look at the large- version of the differential 
equation (|3.3|). Although this equation will acquire additional corrections at multicriticality, as it 



happens at multicriticality for polynomial potentials |17], it will be sufficient to derive the correct 
scaling limit. Assuming also the relation ( 2.33| ) to hold, which will again have to be corrected as in 
p7[|, implies ^(A) ~ A**^^ at A = 0. This leads to the following approximate equation at A f« 



#(A)" 



1 



(4.21) 



A^' 



2 5'^'^^ 



tan'^ 


'KIT' 


^2.-2 / 




I 2 J 





A2 

^2 



+ 



A2 



A2 



. 



It is easy to see that in order to obtain a nontrivial large- A^ limit we have to define 

^ = A^V«A 



(4.22) 



as the appropriate microscopic scaling limit at criticality^. The naive differential equation obtained 
with the rescaling eq. ( [4.22 ) then reads 



tan 



KTT 



^2k-2 ^ (jl)^a« 



a 



. 



(4.23) 



This equation is clearly no longer of Bessel type and thus the corresponding correlation functions 
following from it belong to a new class. However, it cannot be the correct differential equation 
because it is not invariant under the following symmetry on wave functions, 



V'2n (A) 



■ 



(4.24) 



This relation has been shown to hold in |17], independent of the fact whether the potential is polynomial 
or not. However, eq. ( [4.23 ) may hold in an approximate sense as in [17| when in analogy to there the 
condition 1 «; 2^ tan^ [^] ^2«-2 satisfied. 

In order to derive the correct differential equation at multicriticality we would have to proceed 
as in 1 17] and analyze the corrections to eq. ([sT 



from ( [3.lD which are no longer sub dominant at 
criticality. However, because of the nonpolynomial part of the potential this is forbiddingly difficult. 



Following |17] we first would have to derive the auxiliary functions A„(A) and Bn{X), eqs. ( 2.23 ) and 
( p. 24 ) respectively, at finite n. Choosing an explicit example for the potential such as eq. ( [4.1[ ) in the 
beginning of this section one would insert it into the definition ( p. 23 ) and try to evaluate it at finite 
n. The obvious way to do so would be to use the recursion relation eq. ( p. 20 ) and the orthogonality 
of the polynomials. We arrive at 



^„(A) = Nrn{g2 + K dt w{t) 



- A*" 
t2 _ A2 



(4.25) 



^ It generalizes the scaling behavior found in from a different matrix model: ^ = jv''^''*'^'''''^'''*''^'' A, g G N. There, 
multicriticality is reached by adding and tuning an external matrix to a Gaussian model. 
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The integral can only be evaluated at large- using the result eq. ( p. 30 ) and not at finite n as we 
need. Let us point out a second difficulty in the finite- analysis of the string equation ( 2.21| ) which 
determines the recursion coefficients v^i as a functions of the coupling constants. Taking the large-n 
limit and imposing rn±i — > Vn we easily obtain 



n 
N 



(4.26) 



using eq. ( 2.301 ). At n = we obtain the normalization condition eq. (|4.3| ) using 2rN = c, as we 
have mentioned already in the beginning. However, it is unclear how to derive the finite-n version of 
eq. ( 4.26 ) in terms of rn±i, rn±2, etc. It may happen that the string equation contains all recursion 
coefficients from rjv down to tq. 

Let us summarize our findings for the multicritical points at the origin obtained so far. We 
have shown that new classes of multicritical models exist, with p(X) ~ A'^~^ for all real k with 
2m — 1 < K < 2m + 1 and m G N, giving explicit examples. We have shown that at criticality the 
eigenvalues scale with A^^/'^ (eq. ( 4.22[ )) and we have given an approximate differential equation for 
the wave functions, eq. ( [4.231 ). Because of the relations (3.13) and ( 3.14| ) to the chiral ensemble all 
the results obtained in this section immediately translate to the chiral matrix model eq. ( 3.1C| ). In 
the next section we will investigate if the nonpolynomial potentials also lead to new critical exponents 
in the large- A^ scaling limit at the edge of the spectrum. 



5 Recovering usual multicriticality at the spectrum edge 

In this section we investigate the multicritical points at the endpoint c of the spectrum. This is 
particularly interesting since it is known from matrix models with polynomial potentials that here 



m 



the free energy displays critical behavior, !F ^ {g — ^c) "^'^j with the critical exponents 7str 
coinciding with the minimal models One might expect that by adding a nonpolynomial part to 
the potential these exponents could be altered, as has been already suggested in |^]. In particular 
one could hope that because of the presence of a real parameter, k > 1, these exponents could even 
be nonrational, thus representing a one-matrix model representation of such conformal field theories, 
supposedly again a nonunitary one. However, we will find that this is not true in our case. We will 
first expand the density around the endpoint A = c and tune to criticality. It turns out that additional 
zeros occur precisely in the same way as in the polynomial case, with p{\) ~ (c^ — A^) 2 in the simplest 
case. We also calculate the free energy and its derivatives and find the same critical exponent as in 

In order to be most explicit we have chosen to work again with the simplest example, a Gaussian 
plus a real power, eq. (4.1). Since we find the same phenomenon as for the corresponding critical 



polynomial potential we do not expect changes when adding higher order terms and tuning to higher 
criticality. 



We start by expanding the density eq. (O) for the potential ^(A) = \g2X? -|- ^lAI**, eq. (^!T|) , 
around the endpoint of support c. We find that it has an expansion in powers (c^ — A^)2 which 
immediately follows from eq. ( p. 81) for its Freud part: 



1 



/9(A) = ^(52 + 2k(7c"-2/7)(c2_a2)^ + o(^(c2_A2)lj at X^^c^ . (5.1) 
By requiring 

= 52 + 2Kgc''-^H , (5.2) 

we obtain an additional zero at the edge of the spectrum. Together with the normalization condition 
for the density, eq. (|2.12| ) 1 = \g2c'^ + Hgc'^, this fixes both coupling constants 52 and g. We thus 
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have for the critical endpoint Cc in terms of g2- 

Cc = ^ • 5.3 

52 - 2) 

The fact that the density vanishes with the same power as the density of the critical quartic potential, 
does not suffice to conclude that the critical exponents of the corresponding free energies are the same. 
The free energy is given as an integral over the full support and thus may very well behave differently. 
We use here the saddle-point approach to compute the free energy recalling some known facts. First, 
we calculate the free energy in general and then switch to the critical point (of its derivatives) which 
turns out to coincide with eq. ( |5.2D . 

The planar free energy is given by inserting the solution p(A) of the saddle-point equation (|2.4D 
into eq. ( |2.3| ). Using that p{X) satisfies ( |2.4|) (here we put the normalization to unity, a = 1) and 
eliminating the Lagrange multiplier K we obtain 

^ = - / dXp{X)V{X)- [ dXp{X)ln\X\ . (5.4) 



Choosing our specific example eq. ( [4.1D we insert the density in the form of eq. (2.11) into eq. (| 
The following rewriting of integrals, 

dX [ ds = [ ds [ — , (5.5) 

J\^\ Jo J-s C 

turns out to be particularly useful and the remaining integrals are straightforward. We obtain 



2k 2 



C2] (^_2)(K + 1) 2 ('^-2)' 2 4 



+ 4k(k + 2) - + 2) ' ('-'^ 



In order to find the critical point and exponent of the free energy we have to study its derivatives with 
respect to some coupling constant. Since we do not yet impose criticality only one of the coupling 
constants g2 and g is fixed by the normalization condition eq. ( 2.12| ). We choose to keep the Gaussian 
coupling g2 as fixed parameter and vary g and thus c as a function of g. Its derivative can be obtained 
from the normalization condition, eq. ( p. 12 ), keeping g2 fixed: 



g2 + 2KHgc'^ ^ 



52(k-2)(c2-c2) ' (^-^^ 



where we have used again eq. ( 2.12| ) as well as the definition eq. ( ^.31 ). The derivative has a pole 
precisely at the critical value Cc- The first derivative of the free energy thus reads 



dJ' -iHc'' ( 1 (k-2)(k + 1 



,2 „2 



dg g2 + 2KHgc''-^ \ ^ + 2) 32k(k + 2) ^^'^ 

(4k + (2 - k)^/2c2 + 8) , (5.8) 



4k(k + 2) 



where the denominator from eq. ( ^.Tj ) has been canceled due to eq. (|2.12| ). Obviously is analytic 
everywhere for k > 1 and at the critical point for the density eq. (^]^) we can just set c = Cc. The 
second derivative can be calculated in a similar way and we obtain 
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which is again analytic everywhere. For the third derivative we finally obtain 

d^T _ SH^c^'' 1 
dg^ 52 (k - 2) (c2 - c2) 



(5.10) 



which has a singularity at c = Cc (and at k = 2 which we can obviously exclude). This is precisely 
where the spectral density develops an additional zero. We have so far obtained that the transition is 
of third order at Cc given by eq. (|5.2| ) and thus that the corresponding exponent "jstr is negative. In 
order to determine its value we still have to find the exponent e, with which the denominator vanishes 
in 5 - 9c- 

Since we cannot directly solve the transcendental equation ( 2.12| ) for c we will expand it around 
the critical point Cc- We make the ansatz 

= cl + A{g - g,Y + 0{g - g^f^^ , (5.11) 

where the critical coupling gc follows from eqs. ( |2.12| ) and ( |5.2[) to be 

gc ^ . (5.12) 



Inserting the ansatz ( 5.1l| ) into eq. ( 2.12| ) and expanding we find 



= Hg,-iK - 2)A^cr\g - gcf' + Hc^ig - g,) + 0{g - g,f+^ , (5.13) 
which only has a solution for 

E = 

2 



e = - . (5.14) 



We can also determine the proportionality constant 

= ^c^^ , (5.15) 

and we finally obtain 

c2 = cl{l- yiHK^ {g - 5c)^ + 0{g - g^)^) . (5.16) 
Therefore it follows from eq. ( |5.10| ) that 



(g-gc)" (5.17) 



and thus 'jstr = ~^ for the first multicritical point of polynomial potentials [p5|]. 

In summary we have shown that for our prototype example of a Gaussian plus a real power 
exactly the same exponent of and number of zeros of the density at the endpoint c occur, when the 
potential is tuned to criticality. The question is if the same is true for higher critical points, which 
can be obtained by adding higher powers to the potential. From eq. ( |2.8| ) it is clear that the spectral 



density of a polynomial plus a Freud type potential always has an expansion in powers of \J (? — 
and that thus the accumulation of zeros at A = c happens the same way as in the polynomial case. For 
the free energy however there is no other way to calculate its exponent than in the tedious analysis 
as presented above, by taking derivatives and expanding the solution of the normalization condition. 
The method of orthogonal polynomials Q typically used to determine the exponents of T for arbitrary 
criticality m does not easily apply here since we do not have an explicit form of the string equation 
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( ^.21 ) at hand for finite n. But as our example above was arbitrarily chosen (for simphcity) we do 
not expect that new critical exponents occur for higher multicriticality compared to [25|. 

Since we have calculated the free energy for arbitrary couplings we may also ask about the order 
of the phase transition for the critical points of Secti on ^ at the origin. The free energy is analytic 
everywhere away from the critical endpoint Cc eq. ( |5.3| ) which differs from the value it takes for 
criticality at the origin (see eq. ( [4.5D ). However, the fact that J- is thus analytic there does not imply 
a smooth crossover. In principle we would have to compare to the free energy and its derivatives 
computed on the two-cut side at the transition |43]. Our findings in Appendix ^ make this a very 
delicate task. 



6 Conclusions 

We have introduced a new class of one-matrix models which contain nonpolynomial powers in the 
potential. These models do not seem to have an immediate graphical interpretation in terms of 
Feynman diagrams. However, they still share properties of models with polynomial potentials such 
as critical exponents related to Quantum Gravity, while in other aspects they show new features: a 
spectral density vanishing proportional to an arbitrary real power at the origin and a new scaling 
behavior of the eigenvalues. 

We have first determined the macroscopic large- spectral density using saddle-point and orthog- 
onal polynomial techniques, giving explicit examples such as the simplest nonpolynomial extension 
of a Gaussian matrix model. When taking the microscopic scaling limit at the origin we could prove 
that the know massless and massive universality classes for both the chiral and nonchiral ensembles 
were maintained, with the spectral density p{0) providing the universal parameter. 

After tuning the coupling constants of the potential to make p(0) vanish we found a set of new 
classes of multicritical matrix models with p(A = 0) ~ A**"^ for arbitrary real k > 1. Our analysis 
completes the study of such critical models, where previously only multicritical models with a density 
vanishing as an even power were known. Similarly to the findings there for higher critical points 
additional couplings have to be introduced and tuned. Our results fall into classes with 2m — 1 < 
K < 2m + 1 where m G N. For each class of criticality a minimal set of m coupling constants has 
to be tuned and we explicitly gave a realization of such potentials. The scaling behavior for such 
multicritical densities was determined and can be tested on the Lattice in applications to the chiral 
phase transition. 

We would like to emphasize the generality of the multicritical models we have considered. Our 
results also apply to the case when the critical exponents of the underlying physical model differ 
from the mean field values as we can incorporate arbitrary real positive exponents. It would be thus 
very interesting to calculate exactly the correlation functions for such models as a function of the 
real parameter k. Here, we could only derive an approximate differential equation for the asymptotic 
wavef unctions, which enter the kernel of orthogonal polynomials and thus determine all correlation 
functions. An exact differential equation is currently known only for even integer values of n. Another 
open question left for further investigation is the order of the corresponding phase transition in our 
model. While the continuity of k at even integer values suggests a third order transition the infinite 
slope of the first critical density with 1 < k < 2 may still indicate a lower order transition. 

To this aim we have also analyzed the possible critical points of the planar free energy in our model. 
Nonanalytic parts in its derivatives only occurred at the edge of the spectrum. For a nonanalyticity 
to occur at the origin we would have to compare with the derivatives of the free energy coming 
from the two-cut side, which we have not been able to determine so far. When studying criticality 
by the accumulation of extra zeros of the density at the spectrum edge we recovered the critical 
exponents known from polynomial matrix models. A first hope of finding new classes of possibly real 
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valued exponents was not fulfilled. Such models would represent nonrational conformal field theories. 
Obviously more complicated nonpolynomial terms have to be added to the one-matrix model potential 
to achieve such goal. This could be thought of as a generalization of the known two-matrix models 
which do allow for representations of rational theories, with one matrix formally integrated out. 
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A The two-cut solution 

In this Appendix we give the solution for nonpolynomial potentials when the support of the eigenvalue 
density splits into two pieces. Since we restrict ourselves to the case of even potentials we will not have 
to enter the subtleties encountered for such multi-cut solutions where non-perturbative contributions 



may change the meanfield solution [Q. Our derivation follows |^ (see also ||3§]) and since the 
formalism has been presented in detail in Section ^ we shall remain brief here. We will calculate only 
the smoothed macroscopic density as well as the free energy. As has been noted in [p7| ] the higher 
point correlation functions do no longer possess a smooth large- limit. 

The main difference to the single support solution is that the recursion coefficients r„ from eq. 



( ^.20 ) approach two different rather than a single smooth function: 



TN = \{C- {-Ifd) , (A.l) 

where c and d are the endpoints of the support a = [— c, —d\ U [d, c] with < d < c. Under the 
condition eq. ( [A.l| ) with support a it has been shown in [^5| that the wave functions will approach 
the following "smoothed" limits: 

#(A)2 = - ^ ^= e{c^ - x^)e{x^ - d^) , (A.2) 

^ ^ ' ^ ^{c^-X^){X^-d^) 

rNT^ Y^(c2 - A2)(A2 - d^) 



Similar to eq. (|2.3C1|) no restrictions on the potential have been made in [35|. The auxiliary functions 



eqs. ( 2.23 ) and (|2.24| ) can thus easily be obtained: 



AW . IaMX) = fit ' , (A.4) 



Note the alternating sign which is present in both quantities due to eq. ( |A.l| ). We have used already 
the symmetry with respect to the origin by writing intervals over the positive part of a only. Inserting 
these quantities into eq. ( |2.29| ) which is exact for finite we arrive at the following large-A expression 
for the smoothed, macroscopic density: 



p(A) - ^ J^dt . (A.6) 
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The normalization condition for the density can be read off from the string equations [^ ] analogous 
toeq. p:2lD : 



1 




TT 



dt- 



dt- 



V'{t) 



(A.7) 



To make contact to Section B we give the two-cut result for the same simple example, the potential 



V{X) = -g^X^ + g\\r . 



(Ai 



Inserting it into the general result eq. ( [A.6| ) we see that the Gaussian terms cancels and we obtain 



pW = — 

.Id 



J.K-2 



dt 



A 



K-2 



t2- A2 



Xt 



'(c2 - A2)(A2 
(c2 -t2)(t2 _^2) • 



(A.9) 



It is very interesting to compare this result at the transition point, when the two intervals of a merge, 
with the corresponding density eq. (^]^) with a single interval support at criticality. By sending d — > 
it is easy to see that eq. ( A.9 ) leads exactly to A2 times a density of Preud type eq. ( ^.10 ) at shifted 
K ^ K — 2 (see also eq. ( p.36| )). We thus arrive at 



gK 

— tan 
27r 



Kir 



L 2 



A 



+ 



gK{K 



2) He'' 



■k{k — — 3) 



A^ 



A^ 

r.2 



(A.IO) 



where we have again chosen non-integer 1 < k < 3. It is easy to see that the normalization conditions 
eq. (A.7) matches with the normalization condition eq. ( [4.3| ) for the single cut density together with 
the multicriticality condition eq. ( [4.51) . Moreover the two densities eqs. ( |AlO|) and (p]) agree exactly 
at multicriticality ([4.5|), which follows from the general form eq. ( 4.201 ) for m = 1. As an immediate 
consequence on the transition point the two-cut free energy following from eq. (^|J) is identical to eq. 
( ^.6| ) together with eq. ( [4. 51 ) which has been calculated coming from the one-cut side. 

In order to calculate also derivatives of the free energy in the two-cut case we would have to first 
solve the integral in eq. (A.9) in order to be able to determine J- from eq. (| 



These integrals 

are not simple elliptic integrals as encountered in [46|, and the presence of the nonpolynomial terms 
makes such a task fairly difficult. 
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